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ABSTRACT 





A method is described intended for distributed calibration of a probe microscope scanner consisting in a search 
for a net of local calibration coefficients (LCCs) in the process of automatic measurement of a standard surface, 
whereby each point of the movement space of the scanner can be defined by a unique set of scale factors. 
Feature-oriented scanning (FOS) methodology is used to implement the distributed calibration, which permits to 
exclude in situ the negative influence of thermal drift, creep and hysteresis on the obtained results. The sensi- 
tivity of LCCs to errors in determination of position coordinates of surface features forming the local calibration 
structure (LCS) is eliminated by performing multiple repeated measurements followed by building regression 
surfaces. There are no principle restrictions on the number of repeated LCS measurements. Possessing the ca- 
libration database enables correcting in one procedure all the spatial distortions caused by nonlinearity, non- 
orthogonality and spurious crosstalk couplings of the microscope scanner piezomanipulators. To provide high 
precision of spatial measurements in nanometer range, the calibration is carried out using natural standards — 
constants of crystal lattice. The method allows for automatic characterization of crystal surfaces at room tem- 
perature. The method may be used with any kind of scanning probe microscope (SPM). 


FOS 

Counter-scanned images 
Nanometrology 

Surface characterization 
Nanotechnology 





1. Introduction 


Usually, a probe microscope scanner is characterized by three ca- 
libration coefficients K,, Ky, K, representing sensitivities of X, Y, Z 
piezomanipulators, respectively (to take into consideration a possible 
nonorthogonality of X, Y piezomanipulators, an obliquity angle should 
be additionally determined) [1-4]. Because of piezomanipulators’ 
nonlinearity [5-7] and spurious crosstalk couplings, the probe micro- 
scope scanner may be described by the above coefficients only near the 
origin of coordinates, where the influence of the distortion factors is 
insignificant. As moving away from the origin of coordinates, the to- 
pography measurement error would noticeably increase reaching the 
utmost value at the edge of the scanner field (see Section 2.5). 

The problem may be solved by using a distributed calibration, 
which implies determining three local calibration coefficients (LCCs) 
K,, Ky, K, for each point of the scanner movement space, which can be 
thought of as scale factors for axes x, y and z, respectively [8-13]. A 


reference surface used for calibration should consist of elements, called 
hereinafter features, such that the distances between them or their sizes 
are known with a high precision. The corrected coordinate of a point on 
the distorted image of an unknown surface is obtained by summing up 
the LCCs related to the points of the movement trajectory of the scanner 
[11]. 

Both lumped and distributed calibration of the probe microscope 
scanner should be carried out by the data where distortions caused by 
drifts (thermal drifts of instrument components plus creeps of piezo- 
manipulators) are eliminated. Otherwise, the measurements will have 
large errors [2,3,12,14,15]. In the present work, to eliminate the ne- 
gative influence of thermal drift and creep on the distributed calibra- 
tion results, the methods are used of feature-oriented scanning (FOS) 
[14] and of counter-scanning [15]. 

This article completes the series of three articles dedicated to the 
FOS-based distributed calibration. In the first article [11] of the series, 
the suggested calibration procedure has been described in detail and 
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also the working formulae have been derived for correction of the scan 
of an arbitrary surface. The second article [12] of the series is devoted 
to the virtual mode of the distributed calibration. The virtual mode 
implies that, instead of measurement of the real surface of a standard, 
the calibration program performs the “measurement” of an image of the 
standard surface, which has been obtained earlier during conventional 
scanning. The virtual mode is intended for simulation of calibration 
process and validation of analytical solutions found in the first article of 
the series. The virtual mode allows significant acceleration of the de- 
velopment and debugging of the real mode. The real mode being the 
subject of the present investigation implies that the distributed cali- 
bration of the scanner is performed by graphite crystal lattice in real 
time. The real mode permits to acquire the calibration database (CDB) 
that completely characterizes the piezoscanner of the microscope in use 
in the lateral plane. A distributed calibration by the suggested method 
including calibration in the vertical plane is planned to implement at 
the next stage of the method development. 


2. Experimental results 


Distributed calibration of scanner of the probe microscope Solver™ 
P4 (NT-MDT Co., Russia) in the real mode was carried out by atomic 
topography of basal plane (0001) of highly oriented pyrolytic graphite 
(HOPG) monocrystal (Research Institute of Graphite, Russia). HOPG’s 
mosaic spread angle makes 0.8° (graphite density 2.24 g/cm°, purity 
99.999%). In order to minimize thermal deformation of the sample, a 
graphite crystal of small dimensions 2 x 4 x 0.3mm was used. Three 
adjacent carbon atoms (or interstices) forming an equilateral triangle 
ABC were selected as a local calibration structure (LCS) [11]. According 
to the neutron diffraction method, the lattice constant a (actually, the 
length of a side of the ABC triangle) of a natural graphite single crystal 
makes 2.464 + 0.002A (relative error + 0.08%) [16]. 

Graphite surface topography was obtained by scanning tunneling 
microscopy (STM). With microscope Solver, scanning is carried out by 
moving a sample relative to a fixed tip. As the tip, a mechanically cut 
NiCr wire with diameter 0.3mm was used. To protect the microscope 
against floor vibrations, a passive vibration isolation system was em- 
ployed. Moreover, the microscope was housed under a thermoinsula- 
tion hood, which also served as an absorber of external acoustic dis- 
turbances. Neither constant temperature nor constant humidity were 
maintained in the laboratory room where the calibration was im- 
plemented. The typical noise level of the tunneling current in the course 
of the measurements made about 20 pA (peak-to-peak). 

The microscope scanner is a piezotube having four electrodes on the 
outer side and a solid grounded electrode inside. The piezotube was 
fabricated by sintering piezoceramic powder mixture similar to PZT-5A 
(lead zirconate titanate Pb(Zr, Ti)O3, ELPA Co., Russia). The maximal 
unipolar voltage applied to the piezotube X, Y electrodes makes 150 V. 
Between the electrodes of the manipulators of the same name, the tube 
piezoceramic is polarized in mutually opposite directions. X, Y DACs 
have capacity of 16 bits. 

During the raster scanning, the probe movement velocity at the 
retrace sweep was set the same as at the forward trace. Immediately 
before every raster scanning, a “training” of the scanner was carried out 
[15]. The scanner training is a repeated movement along the first line, 
which allows to decrease creep at the beginning of the scan [17]. While 
training, the actual scanning velocity was also determined. 


2.1. Selection of a calibration area on the surface of a standard 


Since the whole scanning field of the microscope used makes about 
2 X 2um and the typical sizes of atomically smooth regions suitable for 
calibration on the surface of the graphite sample did not usually exceed 
450 x 450nm, the distributed calibration was carried out by parts. 
First, a trial scanning of the whole field was conducted, and then plane 
smooth typically 300 x 300nm size regions with no visible defects 
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were selected therein. 

The above scale is set based on the following considerations. First, 
nanodimensional surface defects, whose presence is undesirable, are 
still distinguishable with this scale. Second, with the typical net step of 
5-25 nm, the time of calibration by a region of such sizes will not be too 
long. Third, around the region, a surface area suitable for calibration 
can still be reserved with the sizes that allow the probe to stay within its 
limits even after being shifted by the end of calibration due to drift. 

After the calibration had been completed at each of the selected 
regions, the whole scanner field was shifted to an adjacent area on the 
crystal or the pyrographite crystal was cleft. Periodical graphite 
cleaving not only permits to refresh the surface but also to reduce the 
errors resulted from small variations of the lattice constant and varia- 
tions of angles between the crystallographic directions during the 
subsequent averaging. Those errors arise from imperfections relating to 
the standard itself (see Section 2.6). 

Since during the time of calibration, which may take up from sev- 
eral hours to several days, the standard surface is shifted significantly 
within the field of view of the microscope because of drift, the whole 
scanner field does not require to be shifted frequently. The calibration 
can be continued at either the shifted old surface region or at a new 
region that appeared in the field of view of the instrument. The op- 
erations described above are repeated until most of the scanner field is 
calibrated. 

Fig. 1 shows the areas of distributed calibrations within the scanner 
field. The areas consist of red points designating the LCS positions for 
which LCCs have been found. The sizes of the areas are determined by 
the net step and the number of net nodes. For example, outlined with a 
single frame is the area 73 x 73A calibrated by net of 31 x 32 nodes 
having the step as small as 8 positions (about 2.44 A). Outlined with a 
double frame is the area 396 X 402nm calibrated by net of 37 x 39 
nodes having the step as large as 361 positions (about 11 nm). Some 
areas overlap with each other. The triple frame outlines one of the areas 
where such overlap takes place. In areas like this, it is possible to 
analyze relative variations of the constant lattice of the applied HOPG 
crystal (see Section 2.6). All the calibrations in Fig. 1 have been per- 
formed with the Z manipulator of the scanner extended to the middle of 
its range. 

The presence of some areas not covered by calibration in Fig. 1 is 
not critical as regression (averaging) surfaces will be drawn through the 
obtained LCCs using the least-squares method (see Section 2.4). The 
suggested method allows conducting repeated and additional calibra- 
tions in any area of the scanner field as many times as needed. The 
obtained data are simply added to the common calibration database. 

One placement of the STM head and one tip were sufficient for 
1.5 weeks of practically ceaseless operation. That period is restricted by 
the tip lifetime, which is determined by the slow growth of native oxide 
layer [18]. The tip surface oxidation causes an increase in the tip cur- 
vature radius as well as in the noise level of the tunneling current, 
which leads to the loss of atomic resolution at a certain moment. The 
increase of the curvature radius is caused by tip material volume in- 
crease during oxidation. The increase in noise level is caused by the 
growing number of water molecules adsorbed by hydrophilic oxide 
[19] while working in the air; another reason of the noise is because the 
oxide is mechanically instable [18]. 

Despite of large amount of scanning typical of FOS [14,15], the 
calibration under consideration implies less mechanical wear of the tip 
than usual provided that the bias voltage and the tunneling current 
have the proper values. The fact is that most of the time the tip is held 
within the boundaries of atomically smooth area of the standard, so 
there is a low possibility of the tip collision with highly prominent to- 
pography asperities. Moreover, the tip moves over the surface with a 
small step since the microscope operates near its resolution limit; the tip 
moves at small distances and with low speed in order to prevent large 
creep. All these factors decrease the probability of mechanical inter- 
action between the tip and the surface that could result in notable tip 
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Fig. 1. Areas covered by distributed calibrations. 
The red points “+” designate the LCS positions for 
which LCCs were found. The x and y axes are 
graduated in lateral positions of the scanner. The 
whole scanner field is approximately equal to 
2 X 2um. The calibrations have been conducted at 
the position of the scanner Z manipulator extended 
to the middle of its range. 
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wear (blunting). 

The change in the shape of the tip caused by the above reasons has 
little effect on the accuracy of the calibration, since all atoms/inter- 
stices in the apertures and segments are distorted by the tip in the same 
way, and therefore there is no change in the relative distances between 
atoms/interstices in the LCS. Moreover, unlike calibration by micro- 
objects, the interaction of the complex side surface of the probe with 
the standard surface is completely absent during calibration by atoms of 
a crystal lattice; only the probe apex, which is actually a small cluster of 
several atoms having a well defined form, is involved in the interaction. 

Graphite cleavage between the calibrations was carried out seldom, 
approximately once in 2-3 months. In case of failure to get a stable 
image of good quality by using the constant-current mode, the constant- 
height mode was applied for calibration in the lateral plane. 


2.2. Measurement parameters and modes 


Before the calibration started, the scanner’s Z manipulator was 
moved to the middle of its range by using a coarse approach stage 
whereupon it was held near this position during all the time of cali- 
bration. Since the drift in the vertical plane is tens of times less than the 
drift in the lateral plane, control commands were applied to the Z 
manipulator rarely. 

The STM measurements of the standard surface were conducted 
under ambient conditions at room temperature after the instrument had 
been warmed up [20,21] in active state for hours. The active state 
implies periodical attachment of the microscope probe to the carbon 
atom [11,14] located near the start position (1st net node) of the dis- 
tributed calibration. Beside warming up, the attachment allows for 
monitoring the instrument state — noise level, resolution, and drift ve- 
locity. Thus, by observing the change in drift velocity in this mode, the 
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termination moment of the warming-up stage can be determined. 

The calibration results for each area were put into the common 
calibration database. The calibration was carried out either by LCS 
consisting of carbon atoms or by LCS consisting of interstices [12] de- 
pending on which LCS was closer to the current net node. Such ap- 
proach allows avoiding calibration interruption in case the image of the 
atomic surface becomes accidentally inverted. Moreover, the probe 
path from the current net node into the A feature position becomes 
shorter. 

Fig. 2 shows the initial net (net nodes are designated as “+”), the 
locations of LCSs detected during calibration (designated as “ +”) as well 
as movement trajectory of the probe from feature to feature (solid line). 
The movements within the local rasters while scanning the apertures 
[11] and segments [11] are not shown at the trajectory. In Fig. 2(a), the 
case of calibration with the small step ~ 2.44 A (8 scanner positions) is 
presented and in Fig. 2(b), the case of calibration with the large step 
= 11 nm (361 scanner positions) is presented. The calibration area size 
in Fig. 2(a) is 73 x 73A (in Fig. 1, the area is outlined with a single 
frame); in Fig. 2(b) the size is 396 x 402nm (in Fig. 1, the area is 
outlined with a double frame). In both cases, the measurements were 
taken in the constant-current mode. According to the trajectory given in 
Fig. 2(a), the microscope probe is continually shifted aside from the 
nodes of the initial net due to some drift (cp. with Fig. 2 in Ref. [12]). 
Fig. 2(b) clearly shows that the movements by the net nodes are per- 
formed in full compliance with the trajectory shown in Fig. 1(a) of Ref. 
[11]. 

The waiting time after moving to the next net node was 3s for the 
calibration shown in Fig. 2(a) and 11s for the calibration in Fig. 2(b). 
The number of probe attachments following the pause usually made 
2-7 with the regular scanning and 2—4 with the counter-scanning. As an 
example, three consecutive probe attachments PA1, PA2, PA3 to the 
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Fig. 2. Distributed calibration within area (a) 73 x 73 A (31 x 32 nodes), (b) 396 x 402 nm (37 x 39 nodes) with step (a) 2.44 A (8 positions), (b) 11 nm (361 
positions). The number of detected LCSs (a) 992, (b) 1443 (“” designates LCS locations). Measurement mode: STM, constant-current, Uy, = 50 mV, In = 856 pA. 
During the time (a) 23.8h, (b) 57.3h of calibration, (a) 24174, (b) 48365 attachments were carried out; (a) 2429, (b) 3010 apertures and (a) 22750, (b) 31569 
segments were scanned. Mean lateral drift velocity made (a) 0.147 A/s, (b) 0.149 A/s. The x and y axes are graduated in lateral positions of the scanner. The data 
correspond to the position of scanner Z manipulator moved to the middle of its range. 


feature A located near one of the net nodes are shown in the inset of 
Fig. 2(b). The attachment PA1 is the first attachment performed after 
moving the tip into a net node. The attachment PA2 is carried out by 
the results of scanning-recognition of an aperture with the center in 
point PA1. The attachment PA3 is carried out by the results of scanning- 
recognition of the segment of feature A with the center in point PA2. 
The total number of probe attachments to feature A during the cali- 
bration within the area in Fig. 2(a) was 24174, within the area in 
Fig. 2(b) — 48365. 

In the course of the measurements, the aperture size was 31 X 31 
points (9.2 x 8.8A approx.), the segment size was 21 X 21 points 
(6.1 x 5.9A approx.). The scanning speed in aperture/segment was set 
as 592 A/s for calibration presented in Fig. 2(a) and as 465A/s for 
calibration presented in Fig. 2(b). During calibration, 2429 apertures 
and 22750 segments were scanned in the area shown in Fig. 2(a), 3010 
apertures and 31569 segments — in the area shown in Fig. 2(b). The 
number of averagings of topography height in a raster point while 
scanning the apertures and the segments was 15 and 25 for calibrations 
in Fig. 2(a) and in Fig. 2(b), respectively. Only one of the LCSs detected 
in the aperture, the one nearest to the current net node, was used for 
calibration (LCS ABC, see Fig. 3(a) in Ref. [11]). 

In the inset to Fig. 2(b), PSkp and Skp are preskipping and skipping 
cycles, respectively (for the sake of clarity, the probe movements in these 
cycles are shown slightly idealized). The number of consecutive skipping 
Skp cycles was set to 4. Before the skipping at least one preskipping PSkp 
cycle was performed [14]. The required number of preskipping cycles is 
set automatically. Preskipping is an idle (tuning) skipping required for the 
relative distance measurement process to reach its operating condition (cp. 
trajectories of PSkp and Skp in the inset in Fig. 2(b)). Movement velocity 
between features was set to 22 A/s in both calibrations presented in Fig. 2. 
LCS measurements were performed at moments when the lateral drift 
velocity did not exceed 0.05 A/s by module. In the inset to Fig. 2(b), the 
arrows near the A, B, and C features show the drift direction. Mean lateral 
drift velocity |),| (thermal drift + creep) during calibration by the area in 
Fig. 2(a) made 0.147 A/s, by the area in Fig. 2(b) - 0.149 A/s. The number 
of feature losses fixed during calibration by the area in Fig. 2(a) was 17, by 
the area in Fig. 2(b) — 68. 
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While conducting the measurements, it was found that the larger is 
the speed of probe lateral movement relative to the surface in aperture/ 
segment or at the skipping, the stronger is the creep excited during this 
movement. The strong creep, in its turn, requires setting a longer pause 
and leads to a greater number of attachments inserted automatically 
after the pause. On the other hand, the movement velocity should not 
be set too small since too slow movement leads to low calibration 
productivity and to higher probability of feature loss in case of a strong 
drift. The above scan and skipping velocities have been found experi- 
mentally and are optimal for this particular microscope. 

The larger the step of the initial net, the stronger disturbance is 
introduced into the calibration process at the moment of movement 
from the current node to the next one and the longer time is required 
for this disturbance relaxation. The method retained full functionality 
while conducting distributed calibrations with the step of the initial net 
up to 500 nm. Larger steps were not tested because of lack of suitable 
regions of large enough sizes on the standard surface. 


2.3. Actions of the algorithm when the measurand goes out of the tolerance 
limits 


Because of defects and mostly because of all kinds of instabilities, 
the measured LCS sizes may not fall into the acceptance limits defined 
before the calibration (usually + 8% for linear dimensions and + 5° for 
angular ones). In such situation, the calibration program, in accordance 
with the operator’s choice, can act on the following two scenarios. On 
the first one, after all the possibilities of rescans have been exhausted, 
the algorithm fixes a surface defect in this place, and after that it 
continues calibration in the next net node. In this case, LCCs for some 
scanner positions will be absent in the obtained CDB. 

According to the second scenario, the algorithm moves the probe 
again into the current net node where a new attempt is made to find 
LCCs within the neighborhood of this node. The number of such returns 
to the current node depends on particular measurement conditions at 
the given moment of time. This state of calibration process is re- 
presented as loops on the trajectory. After some time, either the mea- 
surement conditions change so allowing to obtain LCS with parameters 
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Fig. 3. Regression surfaces of 2nd order (upper row) and 1st order (lower row) drawn through (a), (d) LCCs K,., (b), (e) LCCs K,, and (c), (f) local obliquity angles a of 
the obtained CDB. The vertical scale of the surfaces (d)-(f) is intentionally stretched out so to better see the slope of those surfaces. CDB size makes 106683 LCSs. The 
x and y axes are graduated in lateral positions of the scanner. The data correspond to the position of scanner Z manipulator moved to the middle of its range. 


within the limits of the set tolerances, or, due to drift, some other area 
of the standard surface with no defects gets into the current node sur- 
roundings. The results shown in Fig. 1 were mainly obtained with the 
use of the second scenario. 


2.4. Calibration database, constructing regression surfaces 


Figs. 3(a)-(c) show the sought 2nd order regression surfaces built by 
LCCs of the joint CDB of the piezoscanner. To facilitate the analysis, 
Figs. 3(a)-(c) have the same vertical scales as Figs. 3(a)-(c) in Ref. [12], 
where the regression surfaces were obtained in the virtual mode. The 
joint CDB contains information about n = 106683 LCSs. Because of 
disturbing factors and, to a much lesser extent, defects of the standard, 
some LCCs and local obliquity angles have noticeable deviations 
(within the limits of the tolerances set at the calibrations) from the 
expected values. Therefore, before being used the CDB was thinned out 
discarding such values. In order to perform the thinning as efficiently as 
possible, first the expected values should be obtained as close to the 
true values as possible. To do so, the method described in Ref. [15] can 
be employed for example. 

Mean values of LCCs <K,> = 0.304A, <Ky> = 0.290A and ob- 
liquity angle <a> = 0.1° were found near the scanner’s origin of co- 
ordinates (area of analysis 500 x 500 nm) using the regression surfaces. 
Indication that the real mode calibration is done correctly, is the 
proximity of the mean values of LCCs to each other [15] and to the 
initial coefficients A, = 0.306 A, Ay = 0.307 A as well as the proximity 
of the mean obliquity angle to a = 0.4° [15]. Note that the coefficients 
A,, Ay and the obliquity angle a were obtained near the scanner’s origin 
of coordinates during calibration by the image which distortions in- 
duced by thermal drift and creep had been eliminated [15]. 

Because of the appreciable difference between <K,> and <K,> 
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(4.6% instead of 0.3% in Ref. [15]), it could be assumed that for some 
reason the y drift component was not fully corrected in the course of the 
distributed calibration. However, after switching the direction of the 
fast scan in apertures and in segments from x to y, exactly the same 
values were obtained. In order to finally eliminate any doubts con- 
cerning the differences in the coefficients, the sample was rotated by 
90°. Had the drift been completely ruled out in the course of the dis- 
tributed calibration and had the used standard have just insignificant 
uncertainties, the same coefficients would have been obtained over 
again with the new position of the sample. And indeed, after the cali- 
bration, the coefficients have not changed. 

Neither has dependence been detected of the difference in the 
coefficients upon the step sizes of the net. It is also worth noting that the 
difference of the coefficients is not only observed in the neighborhood 
of the origin of coordinates but in other parts of the scanner field and 
for the whole scanner as well. For instance, the mean values calculated 
by the whole CDB are as follows: <K,> = 0.306 A, <Ky> = 0.292 A 
(difference 4.6%), <a> = 0.1°. Thus, the obtained result shows that the 
difference between the coefficients is caused by some systematic error 
which has nothing to do with the standard nor with the method used for 
compensation of drift influence. 

Since the counter-scanned images used for calibration in Ref. [15] 
were obtained more than 8 years ago (the same pyrographite crystal 
was used both then and now), a certain degradation of piezoceramics, 
mechanical or electronic units of the microscope may have probably 
occurred for such a long period of time. On the one hand, taking into 
account that <K,> exactly coincided with A,, we could assume that the 
piezoceramics itself has not changed, since it is hardly possible that 
exactly half of the scanner has aged while the other one has not. On the 
other hand, during the raster scanning, the workload of X manipulator 
is much higher than that of Y manipulator, which could cause the 
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Fig. 4. Static nonlinear distortions of the probe microscope piezoscanner in the 
lateral plane. The arrows show the value and direction of the distortion. The 
length of the vector located in the right top corner of the scanner field makes 
33 nm. The picture corresponds to the position of scanner Z manipulator moved 
to the middle of its range. 


difference in sensitivities due to the difference in aging rates. The 
previous operation of the scanner during many years with the similar 
workload of the manipulators did not lead to any differences in sensi- 
tivity between the X and Y manipulators according to the calibration 
results obtained in Ref. [15]. There were no long breaks in microscope 
operation [22]. More detailed analysis showed that we most likely 
observed a disadjustment of electronics of the piezoscanner’s Y channel 
that apparently happened when transporting the instrument. 

The comparison of the 2nd order regression surfaces in 
Figs. 3(a)-(c) obtained in the real mode with the 2nd order regression 
surfaces in Figs. 3(a)-(c) in Ref. [12] obtained in the virtual mode 
showed that the surfaces obtained in the real mode are close to planes 
(vertical scales in these figures are equal). The detected proximity in- 
dicates that the main causes of the regression surface bending are dis- 
tortions caused by thermal drift and creep, which are excluded by FOS 
methods in the real mode. 

If the fact is established that the microscope scanner can be char- 
acterized by 1st order regression surfaces (see Figs. 3(a)-(c)), then the 
subsequent distributed calibrations of this or similar scanners can be 
substantially simplified. As is known, the position of a plane in space is 
uniquely defined by coordinates of its three points. Thus, by conducting 
the distributed calibration in three positions of the scanner field, lo- 
cated as far as possible from each other, the searched for regression 
planes can be found [23]. To reach the required precision of correction 
of nonlinearities and spurious couplings, the appropriate number of 
repetitive calibrations should be carried out in each of the three posi- 
tions. Since LCCs change very smoothly and the movement range of Z 
manipulator is usually by an order less than the movement ranges of 
X, Y manipulators, the search for LCCs can be carried out for, say, five 
approximately equidistant positions of Z manipulator. Thus, the dis- 
tributed calibration in the space of scanner movements in that case can 
be reduced to measurement of LCCs in only 15 points. 


2.5. Nonlinear piezoscanner distortions and spurious couplings between 
manipulators 


Replacing the found 2nd order regression surfaces with 1st order 
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regression surfaces (planes) and changing the vertical scales so as the 
tilts of the regression surfaces are clearly seen, we obtain Figs. 3(d)-(f). 
The regression surface Ky (see Fig. 3(d)) is tilted towards the axis x, 
which points out a nonlinear response of the scanner X manipulator to 
the driving voltage applied. Besides, the regression surface Ky is also 
slightly tilted towards the axis y indicating small nonlinear spurious 
couplings from other scanner manipulators, in particular, the Y ma- 
nipulator. As is seen from the figure, the response nonlinearity of the X 
manipulator noticeably exceeds the nonlinearity due to spurious cou- 
plings. 

The regression surface Kj, (see Fig. 3(e)) in comparison to the re- 
gression surface Ky has notably less overall tilt (the vertical scales in the 
figures (d) and (e) are the same), i.e., contribution of Y manipulator to 
distortions is notably less than the one of X manipulator. It should be 
noted that the tilt of the regression surface K}, towards the axis x is 
comparable with its tilt towards the axis y. This means that with the Y 
manipulator, the nonlinearities caused by spurious couplings are com- 
parable with piezoceramics nonlinear response. The probable cause of 
the observed differences between Ky and Kj is mechanical and/or 
material unsymmetry of the particular piezoscanner. 

According to Fig. 3(f), the regression surface a” is located practically 
horizontal, therefore the dependence of nonorthogonality between the 
scanner manipulators X and Y on the scanner movements in the lateral 
plane is negligible. For numerical estimates of static nonlinear scanner 
distortions, it is convenient to use the maximum differences 
Age = = 0.007 A, Age = = 0.006 A, and Am* = 0.09" of the regression 
sa ices from the horizontal plane (see Eqs. (13) in Ref. [11]), ie., the 
surface into which the regression surfaces would degenerate in case of 
no distortions. 

In order to better understand the degree and the character of the 
microscope scanner static lateral distortions, a vector field is shown in 
Fig. 4, where the arrows represent local values and local directions of 
the revealed nonlinear distortions. The pattern in Fig. 4 corresponds to 
the case when the entire field 2019.3 x 1907.1 nm is “scanned” ac- 
cording to the raster trajectory of the direct image. Nonlinear raster 
distortions are determined by formulae (8) given in Ref. [11] by using 
the joint CDB. 

As it was expected, the distortions are negligibly small near the 
origin of coordinate system of the scanner (supplying voltages are equal 
to zero) and they increase as the probe approaches the edge 
of the range, where the measurement error reaches A?}*100%/ 


KO, 0) = 2.3%. The size of the longest vector in Fig. 4 mere in the 
right top corner of the field) makes 33.4 nm which is about 1.7% of the 
scanner field side. The obtained error values show that the contribution 
of nonlinearities and spurious couplings to the total raster distortion is 
several times less than the contribution from thermal drift and creep 
which makes 8-11% (see Refs. [12,15]). It can be clearly seen in Fig. 4 
that the vectors originating from points on the y axis are located at 
some angle to this axis pointing out the existence of a small non- 
orthogonality (<a> = 0.1°) [15]. The nonorthogonality is equal to the 
angle at which those vectors are seen from the origin of coordinates [1]. 

To determine the residual lateral error of the distributed calibration 
method, root-mean-square deviations of LCCs are found by formulae 
(14) given in Ref. [11]: og, =0.025A, oz, =0.026A (q = 1.1%. 
Taking into account the number n of LCSs in the obtained CDB, the 
estimated error of regression LCCs could make o = 30z,/Vn = 0.0002 A 
(30). The error o is numerically an order of magnitude less than the 
amplitude of thermal motions of carbon atoms in graphite lattice sites 
at room temperature [24]. 


2.6. Errors caused by defects of the standard surface. Automatic 
characterization of a crystal surface 


Since a crystal lattice is used for LCC determination, an additional 
error is introduced in the LCCs because of a variation of the lattice 
constant caused by random imperfections of the standard. Another 
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component of the error o is noises of the registration and positioning 
systems of the STM. With CDB as large as n = 90000 LCSs, the error 
0, = 0a/<K,> = 0.002 A (30) of STM-measurement of graphite lattice 
constant a, caused by the error o of scanner calibration, becomes 
comparable with the measurement error of graphite lattice constant by 
X-ray diffraction or neutron diffraction methods [16]. Unlike the last 
two methods, where the lattice constant is mainly measured within the 
crystal volume, the suggested method enables accurate measurement of 
the lattice constant exclusively on the surface. 

In different locations of the scanner field, the LCC root-mean-square 
deviations are numerically close to each other and to the LCC root- 
mean-square deviations ox,, ox, determined for the whole field of the 
scanner. The relations detected point out that in spite of localization of 
the defect in graphite, random changes of the graphite lattice constant 
are most likely small and the main part of the errors ox,, ox, relates to 
the noise components of the registration and the positioning systems of 
the used microscope. 

Graphite lattice can have small variations of sizes in different places 
of the crystal surface or in the same place before and after cleavage. 
Those variations were numerically estimated by comparing the LCCs 
K1,, K2, (K1,, K2,) from CDB that corresponded to the same area of the 
scanner field (in Fig. 1, one of such areas is outlined with a triple frame) 
but had been taken either in different places 1 and 2 of the standard or 
in the same place but on different cleavages 1 and 2. This approach 
eliminates the influence of the scanner nonlinearities and allows the 
relative change 6, of the lattice constant to be estimated by the relative 
change 0g, = 2|K1,-K2,|100%/(K1, + K2,) (or dx,) of the LCC. For the 
lattice constant located in parallel to X manipulator, 6, = 6x,; for the 
lattice constant located in parallel to Y manipulator, 6g = 6x,. 

A preliminary analysis has shown that the relative change in con- 
stant lattice 6, on the surface of the pyrographite crystal in use makes 
Ox, ~ Ox, = 0.8% which is 5 times higher than in the volume of natural 
single crystal graphite [16]. Thus, the method suggested enables not 
only for scanner calibration but also for estimating errors of the 
standard itself. Matching og, = 0.025 A with the difference 
K1,-K2, = 0.003 A confirms the earlier made assumption that random 
changes in the pyrographite lattice constant are insignificant in com- 
parison with the noise components of the registration and the posi- 
tioning systems. Because of the importance of the issue of invariability 
of HOPG crystal lattice constant, a special study is to be conducted in 
order to cover the issue in detail. 

By saving the relative coordinates of B and C features (skippings 
A<B, A<C) found during the calibration, it is easy to obtain a dis- 
tribution of crystal lattice constants and angles between crystal- 
lographic directions over the sample surface [23]. Thus, beside cali- 
bration, the developed method can also be used for characterization of 
crystal surfaces, localization of surface/subsurface defects, mechanical 
stresses, studying superstructures (superlattices) [25], surface re- 
constructions, adsorption/desorption processes, etc. In particular, by a 
small periodical modulation of sizes of unit cells, low contrast super- 
structures [25] could be discovered, which is hard or impossible by 
using other methods. 

When implementing the surface characterization, instead of rigid 
binding of the net to the absolute coordinates of the scanner; it is ne- 
cessary to calculate each next absolute net position relative to the po- 
sition of the feature A determined in the end of the last skipping cycle 
A <C (see inset in Fig. 2(b)). In this case, the microscope probe will 
drift along with the surface exactly like in FOS [14,25,26] and it will 
never leave the limits of the crystal area under characterization. 


3. Discussion 


At present, microscopes equipped with a closed-loop positional 
system are unable to measure surface topography with high resolution 
because of noises of the linear position sensors. To operate microscope 
near or at the utmost resolution, the closed-loop positional system has 
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to be turned off. As a result, the measured topography is distorted by 
drifts, nonlinearities, and spurious couplings of the scanner manip- 
ulators. 

To correct the errors caused by the manipulators’ nonlinearities and 
spurious couplings, with the closed-loop positional system turned on, a 
sort of distributed scanner calibration is carried out by the standards 
which do not require the utmost resolution of the microscope [9,10]. 
The data obtained during such distributed calibration, which describe 
the nonlinear scanner behaviour, are used later to correct the topo- 
graphy measured with a high resolution. 

The fundamental distinction of the suggested distributed calibration 
method is that the random velocity fluctuations of the drift thermal 
component have no influence on the obtained LCCs. Moreover, the 
error of LCC determination can be substantially reduced by performing 
repeated scanner calibrations the number of which has no principle 
restrictions. Unlike the closed-loop positional system, the developed 
calibration approach does not introduce any additional noise which 
permits calibrating the scanner by natural references such as constants 
of crystal lattices, using separate atoms and interstitials as features. 

Thermal drift of a microscope can be decreased by means of mi- 
croscope temperature stabilization in a special thermally-stabilized box 
or room. The drawback of temperature stabilization is that equalizing 
the temperature of different parts of the instrument takes a rather long 
time and cannot be done completely because of the presence of varying 
heat sources inside the microscope. Temperature stabilization of the 
microscope at a low temperature (2-20 K) makes it possible to further 
reduce the thermal drift rate and achieve high mechanical stability of 
the tunnel gap [27]. However, cooling to a low temperature would 
greatly complicate the experiment and increase its conduction time. 
Moreover, the cryogenic probe microscope is a much less universal 
instrument as compared with regular microscopes, which substantially 
restricts its application area. 

Due to the built-in mechanism of in situ compensation of drift in- 
fluence [11], the distributed calibration and the crystal surface char- 
acterization can be performed by the suggested method in nanometer 
and subnanometer scales at room temperature without special tem- 
perature stabilization of the instrument. Not high mechanical stability 
of the tunneling gap at the room temperature is compensated by the 
possibility of the calibration method to perform multiple repeated 
measurements. 


4. Conclusion 


The use of the atomic surface of a crystal as a standard surface al- 
lows obtaining LCC distribution of a very high density. However, at 
present it is impossible to cover the whole scanner’s movement range 
with such a dense net, especially at several levels, because of low 
performance of the existing microscopes. The results obtained in the 
course of real distributed calibration show that the dense net should not 
necessarily be used. For example, to acquire an adequate CDB, it is 
enough to perform calibration only in three points of the scanner field 
for each of 3-5 equidistant positions of the Z manipulator. 

The key points of the research can be summarized as follows: 


(1) Operability of the distributed calibration method built on the base 
of feature-oriented approach has been confirmed. 

(2) Feasibility of accurate calibration of the probe microscope scanner 
by a crystal lattice of highly oriented pyrolytic graphite at room 
temperature has been proved. 
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